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We report preliminary results from a study of the decay B" -^ D'^tv'^tt^ using a data sample of 
470.9 ± 2.8 million BB events collected with the BABAR detector at the T{4S) resonance. Using the 
Dalitz-plot analysis technique, we find contributions from the intermediate resonances _D2(2460)~, 
Do(2400)~, p{770)'^ and /2(1270) as well as a -k'^tt~ S-wave term, a D^tv~ nonresonant S-wave term 
and a virtual D" (2010) ~ amplitude. We measure the branching fractions of the contributing decays. 

INTRODUCTION 

The study of the Dalitz plot 'l| of B^ — >■ D'^tt+tt^ decays is motivated by several factors. The branching fractions 
of B — )■ D** transitions^ are of interest to help address a conflict betvifeen theoretical predictions [2| and experimental 
results [3[y] in semileptonic B — > D**li> decays. The ZJ^'tt+tt^ final state allows relatively clean studies of the J^ = 0+ 
and 2+ D** states, since the 1+ mesons cannot decay to Dtt. Measurements of these decays test theoretical models 
including quark models [5|, QCD sum rules [6|48| and lattice QCD [9j. Similarly, measurement of the branching 
fraction of the B'^ — > D'^p'^ decay will help to test the dynamics of "color-suppression" in B decays (related to the 
fact that the color quantum numbers of the quarks produced from the virtual W boson must match that of the 



^ D** mesons are P-wave excitations of states containing one charmed and one light (u, d) quark. 



spectator quark in order for a p° meson to be formed) [10l4l3|. Moreover, using isospin symmetry to relate the decay 
amplitudes of i3° — >■ D^p^, B^ — > D~ p'^ and 5+ — )■ D^p'^, it is possible to study effects of final state interactions in 
these decays III llJ] . 



Another motivation is that the S" — >■ D^p^ decay can be used to measure sin2/3, where j3 is the CKM unitarity 
triangle angle la [iGJ , if the D^ meson is reconstructed in a CP eigenstate. The measurement of this angle in the 



b — >■ cud quark-level transition is theoretically cleaner than the commonly used b — ^ ccs decays (such as B'^ 
J/ipKg) [171, |l8| and comparisons of the values measured in different quark- level transitions can be used to search for 
the influence of physics beyond the Standard Model [19J. The time-dependent analysis of the B^ — >■ D^tt+tt" Dalitz 
plot not only allows a proper handling of effects due to interference between broad resonances, but also enables an 
improved measurement of /3 since terms proportional to cos2/3 as well as am2/3 can be measured 20|, l21| . For such an 



analysis, it is necessary to have a good understanding of the population of the B^ — >■ D'^tt+tt" Dalitz plot. This can 
be best studied in the Z?" -^ K^n^ decay, which is the subject of this study. 

The B^ — > £)°7r+7r~ decay has been previously studied by Belle [22| and the related B^ — )• D^n^n^ decay has 
been studied by both BaBAR [23| and Belle [24)]. In this paper we present preliminary results from the first study 
of the B^ — )■ D^-K~^Tr~ decay by BaBAR. The data used in the analysis, collected with the BaBAR detector [25| at 
the PEP-II asymmetric energy e+e" collider at SLAG, consist of an integrated luminosity of 429 fb~^ recorded at 
the T{4:S) resonance ("on-peak") and 45 fb^^ collected 40 MeV below the resonance ("off-peak"). The on-peak data 
sample contains the whole BaBAR dataset of 470.9 ± 2.8 million BB events. 

SELECTION 

We reconstruct _B" -^ D^tt+tt" candidates (the inclusion of charge conjugate reactions is implied throughout this 
paper) by combining a D" candidate with two oppositely charged pion candidates. The charged pion candidates 
are required to satisfy particle identification requirements that have efficiency above 97 % and kaon misidentification 
probability below 20%. We reconstruct D'^ mesons in the decay channel K^tt~. For the Z?" daughters, the charged 
kaon candidates are required to satisfy particle identification requirements that have efficiency above 97 % and pion 
misidentification probability below 15 %, while the charged pion candidates are required to pass slightly looser criteria 
than those for the bachelor pions. The D candidates are required to have an invariant mass within 15 MeV/ c^ of the 
nominal D^ mass [26|; this requirement is 85% efficient for signal Monte Carlo (MC) events. 

The D^ candidate and the two bachelor pion candidates are required to originate from a common vertex. Signal 
events are distinguished from background using two almost uncorrelated kinematic variables: the difference AE 
between the CM energy of the B candidate and -\/s/2, and the beam-energy-substituted mass mES = \/s/4 — p'^, 
where ^/s is the total CM energy and ps is the momentum of the candidate B meson in the CM frame. We apply 
preselection criteria of -0.075 GeV < AE < 0.075 GeV and 5.272 GeV/c^ < mEs < 5.286 GeV/c^; these requirements 
are 86 % efficient for signal MC events. We make further use of these kinematic variables to discriminate signal from 
background in the fit described below. We exclude candidates consistent with the abundant B'^ — > Z?*(2010)~7r+ 
decay by rejecting events which contain a candidate with D^Tr~ invariant mass below 2.02 GcV/c^ (to maintain the 
symmetry of the Dalitz plot, we also remove the region with D'^tt^ invariant mass below the same value). These 
events are used as a control sample to monitor differences between data and MC. 

To suppress the background contribution from continuum e+e~ -^ qq {q = u, d, s, c) events, we construct a neural 
network (NN) discriminant that combines four variables commonly used to separate jet-like qq events from the more 
spherical BB events. These are: the 0*^ order momentum-weighted monomial moment,^ Lq; the ratio of the 2"*^ order 
momentum- weighted monomial (L2) to that of O"^ order (Lq), L2/L0; the absolute value of the cosine of the angle 
between the B direction and the beam (z) axis, |cos0Bmom|; and the absolute value of the cosine of the angle between 
the B thrust axis and the beam (z) axis, |cos0Bthr|. AH these variables are evaluated in the e~^e~ center-of-mass 
frame. We apply a requirement on the NN output that retains approximately 88 % of the signal and rejects ~ 52 % 
of the continuum background. Most of the remaining background originates from B decays, and is discussed below. 

After applying all selection criteria, we retain 26334 events with candidate 5° — >■ ZJ^tt+tt" decays. Around 20% of 
these events have multiple candidates. When an event has multiple candidates we retain the candidate with the best 
geometrical i?- vertex probability. 



The momentum- weighted monomial moments are defined Lj = JZiPi \cos 6i\-' , where 9j is the angle of the track or neutral cluster i 
with respect to the signal B thrust axis, pi is its momentum, and the sum excludes the daughters of the B candidate. 



The efRciency for signal events to pass all the selection criteria is determined as a function of position in the Dalitz 
plot (DP). Using a Monte Carlo (MC) simulation in which events uniformly populate the phase-space, we obtain an 
average efRciency of approximately 35%. The efRciency is shown as a function of phase-space in Figure [TJ both in 



terms of the conventional DP (for which we choose axes m 
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and 
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Dalitz plot" (SDP). The latter is described by the variables M and 0, 
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where mjr+Tr- is the invariant mass of the two pions, 



tobo - m^o and m^^^ 



and Ott+tt- is the angle between the 1?° and the 7r+ in the tt 
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_ — 2to^ are the kinematical 
rest frame. While the conventional 
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DP representation provides a useful visual representation of the physics of the signal decay, the SDP allows closer 
scrutiny of the most densely populated regions of the phase-space, and hence is appropriate for studies of background 
distributions, for example. 
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FIG. 1: Variation of the signal reconstruction efficiency over the phase-space. In these plots the D*(2010) veto is not applied. 



BACKGROUNDS 

In addition to the background from continuum processes, we expect backgrounds from other BB decays. These 
are studied using large MC samples in which the B mesons decay generically according to our current knowledge 
of their branching fractions. We classify backgrounds from BB decays in six categories based on their /S.E and 
Toes distributions as determined from MC samples. The different BB background categories also have different DP 
distributions. 

Table H] lists the expected number of events and the dominant contributing mode for each category. Background 
categories 1 and 2 have four track final states, and peak in both IS.E and toes ~ category 1 has signal-like peaks in 
both, while category 2 has a IS.E peak shifted to positive values due to pion— ^kaon misidentification. The decay modes 
that contribute to these categories do not contain real D mesons, with the exception of D^K^^, which contributes to 
category 1. Background categories 3-6, which are dominant, do contain real D mesons. Category 3 peaks in toes 
and has a /S.E distribution that is shifted to negative values due to kaon— ^pion misidentification. Categories 4-6 do 
not peak strongly in either /S.E or toes- Category 4 has a broad toes distribution and a slight peak in A£', and 
includes background from £)*(2010)~7r+ events that escape the veto due to misreconstruction. Category 5 has a 
broad ttjes distribution (similar to category 4) and an approximately linear /S.E shape. Category 6 has combinatorial 
distributions for both tties and /S.E. The continuum background shape is combinatorial and does not peak strongly 
in either IS.E or tties . A summary of the backgrounds in given in Table ID 



MAXIMUM LIKELIHOOD FIT 



We perform an extended unbinned maximum likelihood fit using the variables /S.E, toes and the DP co-ordinates 
in order to determine the signal yield and the properties of the Dalitz plot. The complete likelihood function is given 



TABLE I: Summary of backgrounds. For each category the dominant contributing mode and the expected number of events 
after all selection requirements are applied to the data are given. 



Category 


Dominant contribution 


Total # Expected 


BB 1 


J/iPK+n' 


444 ± 24 


BB 2 


a^TT^F 


32 ±7 


BB 3 


D°K+7v- 


240 ± 18 


BB 4 


D°p+ 


7415 ± 101 


BB 5 


D*°7V + 


1475 ± 44 


BB 6 


Combinatoric 


7336 ± 99 


qq 




5352 ± 226 



by: 



exp 



N, 



-E^^ n 



J2NkVk{ml\mV,mEs\AE') 



(2) 



where Nk is the event yield for species fc, the index i runs over the N,, events in the data sample and Vk is the 
probability density function (PDF) for species k, which consists of a product of the DP, toes and AE PDFs. The 
different species k are signal, qq background and six BB background categories. The function — ln>C is minimized to 
obtain the preferred values of the free parameters of the fit. 

For each of the BB background categories, the AE, toes and DP PDFs are described with histograms obtained using 
MC. For qq background, the AE and ttt-es PDFs are a l'^*-order polynomial and an ARGUS function \^7}, respectively. 
The parameters of the ARGUS function are fixed to values determined using off-peak data, while the slope of the qq 
AE PDF is a free parameter of the fit. The continuum background DP PDF is modelled with a histogram obtained 
from data in a sideband region of toes , after subtraction of the (MC-based) expected contribution from BB decays in 
this region. We have verified the consistency of our background PDFs in off-peak data, in background MC samples, 
and in on-peak data sidebands. All histograms used in the fit are in the square Dalitz plot format. 

The signal component is composed of two parts which are distinguished by whether or not the kinematics of the 
daughter particles are well reconstructed. We refer to the well reconstructed events as "correctly reconstructed" (CR) 
and the misreconstructed events as "self-cross-feed" (SCF). The fraction of SCF events as a function of DP position 



/scf(w 



+ '' 



J is determined from MC, and is shown in Figure [51 Its value is typically below 10 % but is larger in the 



corners of the Dalitz plot where one of the pions has low momentum. 
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FIG. 2: Fraction of misreconstructed events as a function of the phase-space. In these plots the D* (2010) veto is not applied. 



Both CR and SCF events have the same underlying physics PDF, but due to misreconstruction SCF events have 
reconstructed DP positions that differ from their true values. This smearing is implemented by convoluting the PDF 



with a resolution function RgQp(m~^_, 



\fh\, 



rh_) that gives the probability that an event with true DP position 



(m^, m._) is reconstructed at (to,^, to._), and is described by a histogram in the square Dalitz plot co-ordinates that is 
itself a function of position in the phase-space. For correctly reconstructed events, DP resolution effects are negligible. 



The signal Dalitz plot PDF is thus written as 



DP 



-pphys ^fjj2^ ^ m^_)e{'m\_ , m^ )/scf ("t-+ , '>n^-)Rsc¥ (™+ , ™- ; ™+7 "t-- ) d{fk\)d{rh'^_ ) 



where 7'P'^^''(r7i^, m?_) is the underlying physics PDF (discussed below), e{'m\^rn?_) is the efficiency (Figure [T]), and 
/scf('ti+,''t^?-) is the SCF fraction (Figure [2|). The integral is over the Dalitz plot. The normalization factor M 
ensures that Psig (m^ , to?_ ) gives unity when integrated over the phase-space. 

The CR and SCF signal events have different distributions in Ai? and mES- For mES, both CR and SCF PDFs 
are described by double Gaussian functions where the widths of the two Gaussians are constrained to be the same. 
For IS.E, the CR PDF is again a double Gaussian function (in which the two Gaussians have different widths) while 
the SCF PDF is represented by an asymmetric Gaussian with power-law tails. The two Gaussian widths of the Ai? 
PDF for the CR component are given by linear functions of (ni^™) = min(TO^, m?L) to account for the momentum 
dependence of the resolution across the DP. All SCF PDF parameters are fixed to values determined from MC, while 
CR PDF parameters are floated in the fit where possible. The CR PDF parameters that cannot be determined from 
the fit are determined from MC. Data/MC correction factors determined from the _D*(2010)^7r+ control sample are 
applied to all such parameters, except for the slopes of the dependence of the /S.E widths on {mJf,™y 



\2 



We determine a nominal signal DP model using information from previous studies of B^ — > D'^tt'^tt 22| and 



B^ -^ D 7r+7r+ [23|, |2j|, and the change in the fit likelihood value observed when omitting or adding resonances. 



We use the isobar model |28l - l3(j | , which models the total amplitude as resulting from a sum of amplitudes from the 
individual decay channels: 

N 

rP''y%ml,ml)^\A{ml,m'i)\'^ where yt(TO2_^ m?_) = ^ Cj-F^ (m^., m^) , (4) 

where Fj (rn^ , rn^ ) are the dynamical amplitudes and Cj are complex coefficients describing the relative magnitude 
and phase of the different decay channels. All the weak phase dependence is contained in the Cj coefficients, which 
we express in terms of their real and imaginary parts: Cj — Xj + iyj, so Fj{rri^, m?_) contains kinematics and strong 
dynamics only. We treat the D" -^ K^tt^ decay as flavour-specific and neglect contributions from the doubly- 
Cabibbo-suppressed D" — >■ K'^tt" decay. We assume direct CP violation is negligible and hence use the same model 
for B^ — > Z3°7r"*"7r~ and its conjugate decay. We also neglect possible contributions from b —>■ u mediated, and hence 
highly suppressed, transitions {e.g. B° — >• D^(2460)^tt^). 

In the Dn spectrum previous studies [22h24{ have observed contributions from 132(2460) and Dq{2A00), as well as 
the effect of a virtual D* (13,* (2010)) amplitude. The latter amplitude is described as virtual since although the region 
around the narrow 13* (2010) pole is vetoed, off-shell production can contribute to the amplitude - the effect is similar 
to a nonresonant P-wave term. We find that an additional nonresonant (S-wave) Dn contribution is necessary to fit 
the data; we describe the nonresonant (NR) term using an empirical shape, first introduced in Ref. [3l| . proportional 

2 

to e~*""^-, where the shape parameter is determined from the data to be a = 0.60 ± 0.15 (statistical uncertainty 
only). In the tt+tt" spectrum previous studies [22j have observed contributions from p(770)° and /2(1270). We find it 
is necessary to include S-wave terms and hence include a contribution using the K-matrix formalism [32h34| , described 
in more detail in the Appendix. To our knowledge, this is the first use of the K-matrix formalism in B meson decays. 
All other resonances are described using relativistic Breit-Wigner (RBW) shapes, with Blatt-Weisskopf barrier form 
factors [35| and angular distributions given in the Zemach tensor formalism 36|, [37 1 . The Dalitz plot formalism used 



in this analysis is the same as that described in more detail in several previous publications 38144 ij. The masses 
and widths of all resonances are constrained to world-average values [26|, while K-matrix parameters are fixed to the 
values tabulated in the Appendix. 

In total there are 43 free parameters of the fit. These are the yields of signal, qq and the 6 BB background 
categories; the real and imaginary parts of 5 intermediate contributions to the signal DP model (not counting those of 
Z32(2460)~7r+ which are fixed as a reference); the real and imaginary parts of 10 complex coefficients in the production 
vector of the K-matrix parametrization of the tt+tt^ S-wave; 2 parameters each of the CR signal AE and toes PDFs 
and the slope of the continuum AE distribution. 



RESULTS 

The fit returns 5098 ± 102 signal events. For tliis and all other quantities the statistical uncertainties are calculated 
from an MC study where the events are generated from the PDFs and the PDF parameters are the central values 
from the fit to data. Yields of the various background categories are broadly in line with expectation, although there 
appears to be some cross-feed between BB categories. Projections of the fit result onto toes and AE are shown in 
Figure [31 while projections onto each of the two- particle invariant masses are shown in Figure H] and projections onto 
the cosines of the helicity angles, defined as the direction of one of the two daughters of the resonance relative to the 
direction of the third particle in the rest frame of the resonance, are shown in Figure [SJ The signal distribution across 
the phase-space, in both conventional and square Dalitz plot co-ordinates, calculated using the sVlot technique J42|. 
is shown in Figure HI Structures due to the D2(2460)^, p(770)° and /2(1270) resonances are clearly visible. 

Figures S] and [5] show that our DP model gives an excellent representation of the data in most regions of the Dalitz 
plot. The only region where discrepancies between the data and the fit result are apparent is at low values of 77i_, 
where a sharp rise near threshold is observed. This structure also appears as a reflection in the ?n+ and cos 0_ 
distributions. We discuss this further when we consider model uncertainties, below. 
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FIG. 3: TTiES and A_E distributions. The points with error bars show the data, the red dotted lines show the continuum 
background, the green dashed lines show the total background, the black dot-dashed lines show the signal, and the blue solid 
lines show the total fit result. 

We calculate the fit fractions and interference fit fractions, shown as a matrix in Table HIl The fit fractions are the 
elements along the diagonal, and are given by 



FF, 



j'^^\cjFj{m\,m?_)\ d{'m\)d{m?_) 



(5) 



while the interference fraction are the off-diagonal elements and are given by 

/j3p2Re [c,c*F,{m\,m'l)F*{mX,mV)\ d{mX)d{ml_] 



FF,, 



/dp J2iCjFj{ml,ml) d{Tnl)d{ml) 



(6) 



for i < j only. Note that, with this definition, FFjj — IFF j. These give a convention independent representation of 
the population of the DP. Although the sum of fit fractions can be greater than unity - in this case it is (148 ± 5) % 
(statistical uncertainty only) - the sum including interference fit fractions must be identically equal to one. The 
largest interference effect is between I?Q(2400)~7r+ and the D-k nonresonant amplitude. 

In Table Hm we give results for the branching fractions. The inclusive i?" — ?> ZJ^tt+tt^ branching fraction is calculated 
by dividing the signal yield by the average efficiency determined from the nominal model, by the number of BB pairs 
in the data sample, and by the branching fraction for the D decay {B{W -^ K~^tt~) = (3.91 ± 0.05) x 10~^ [26i]). 
The average efficiency is found to be 30.6 % and is further corrected for the measured data/MC differences (discussed 
under systematic uncertainties below). Our result compares well to that of Belle: B{B'^ — >■ D'^tt+tt") = (8.4 ± 0.4 ± 
0.8) X 10~* [24I. The product branching fractions for the contributing decay modes are obtained by multiplying the 
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FIG. 4: Projections of the fit results onto the invariant masses. The top row shows the projection of all the data; the bottom 
row has the signal component enhanced by additional requirements (5.276 GeV/c^ < rriES < 5.282 GeV/c^ and |A_E| < 20MeV). 
The points with error bars show the data, the red dotted lines show the continuum background, the green dashed lines show 
the total background, the black dot-dashed lines show the signal, and the blue solid lines show the total fit result. 



TABLE II: Matrix of fit fractions and interference fractions (central values only without uncertainties). 





-D: 


(2460) -7r+ 


^c 


(2400) -7r+ 


p(770)"D" 


/2(1270)D" 


D, 


(2010)-7r+ 


Dtt NR 


K matrix 


D2(2460)-7r+ 




0.2047 


















DS(2400)-7r+ 




0.0000 




0.2481 














p(770)°Z?« 




-0.0133 




0.0264 


0.3343 












/2 (1270)1)'' 




-0.0130 




0.0223 


0.0000 


0.0983 










D;(2010)"7r+ 




0.0000 




-0.0001 


-0.0565 


-0.0347 




0.1579 






Dn NR 




0.0000 




-0.2471 


-0.0246 


-0.0458 




0.0001 


0.1844 




K matrix 




0.0019 




-0.0672 


0.0000 


-0.0003 




-0.0016 


-0.0303 


0.2559 



inclusive branching fraction by the relevant fit fraction. Where possible, these have also been corrected for subdecay 
branching fractions (e(p(770)° -> tt+tt-) = (98.9 ± 0.16)%, S(/2(1270) -^ tt+tt^) = (84.8+?;^)% [1^). We are not 
able to perform such a calculation for 1)2(2460)^71+ since, although decay modes other than Dtt have been seen, 
the relative branching fractions are not knovifn. The Z)o(2400) has only been observed to decay into Dtt, but it may 
be presumptuous to conclude that its branching fraction is 100%. Our results for D2(2460)^7r+, p(770)°I?*' and 
/2(1270)Z)° are consistent with those of Belle, while we see a somewhat larger branching fraction for Z?*(2010)^7r+ 



and a much larger branching fraction for _Dq(2400) 7r+ (Belle measures B{B^ — >■ i)g(2400)" 



X 6(1?;; (2400)- 



D" 



(0.60 ± 0.13 ± 0.15 ± 0.22) x lO^"*). 



SYSTEMATIC UNCERTAINTIES 



We consider the following systematic effects on the values of the fit fractions. 

• Fixed shapes of the efficiency, qq and BB Dalitz-plot histograms: 

The contents of all bins of square Dalitz plot histograms used to describe these shapes are fluctuated in accor- 
dance with the uncertainties. This procedure is repeated many times and the RMS of the distribution of the 
change in the fit results is taken as the associated systematic uncertainty. 

• Fixed TTiEs and IS.E PDF parameters (or histograms): 

We vary any fixed parameters in the PDF descriptions by their uncertainties, taking correlations into account. 
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FIG. 5: Projections onto cosines of helicity angles. The signal component has been enhanced in all plots by additional cuts 
(5.276 GeV/c^ < tties < 5.282 GeV/c^ and |A£;| < 20MeV). The top row shows the projection onto the cosine of the tt+tt" 
helicity angle in the regions (left) around the p(770)'^, and (right) around the /2(1270). The bottom row shows the projection 
onto the cosine of the Dtt helicity angle in the regions (left) below and (right) around the D2(2460). The points with error 
bars show the data, the red dotted lines show the continuum background, the green dashed lines show the total background, 
the black dot-dashed lines show the signal, and the blue solid lines show the total fit result. 
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FIG. 6: sVlots of the signal distribution in the (left) Dalitz plot and (right) square Dalitz plot. 



The variation in the fit results is taken as the systematic uncertainty. For most parameters, their values and 
uncertainties are determined from data control samples. An exception is the self-cross-feed fraction which is 
obtained from Monte Carlo. To conservatively allow for possible data/MC differences in the behaviour of the 
SCF component, we apply a Dalitz- plot independent scale factor that alternately increases and decreases the 
SCF fraction by a factor of two, and take the larger difference compared to the nominal result as the uncertainty. 
The contents of the histograms used to describe the BB background m-Es and AE PDFs are varied using the 
same prescription as described above. 

• Fit bias: 
We generate large ensembles of pseudo-experiments, containing fully simulated signal events, using the param- 
eters returned by the fit to data. From the distribution of results of these ensembles, we evaluate biases on the 
fit parameters. All biases are found to be small compared to the statistical uncertainties. We assign systematic 
uncertainties of the sum in quadrature of half the bias and its uncertainty. 
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TABLE III: Branching fraction 
third is due to the DaUtz-plot 
column gives the product of the 
intermediate resonance decay to 



results from the fit to data. The first uncertainty is statistical, the second is systematic, the 
model, and the fourth (where present) is due to secondary branching fractions. The third 

branching fraction of the B decay to the mode listed in the leftmost column with that of the 

the final state particles. 



Resonance 



Inclusive B 



TT 



fcr: 



Fit Fraction 
(%) 



BiB" -^ Mode) 
xBiR^hh) (10"'') 



B{B" -^ Mode) 
(10-*) 



D"tt 



8.81±0.18±0.76±0.78±0.11 



D2(2460)-7r+ 


20.5±0.9±1.3±3.7 


DS(2400)-7r+ 


24.8±2.5±3.0±12.9 


p(770)"D'' 


33.4±2.0±5.2±10.0 


/2(1270)L>° 


9.8±1.1±1.6±3.4 


D:(2010)-7r+ 


15.8±0.9±1.2±3.7 


Dtt nonresonant 


18.4±2.3±4.3±13.6 


K matrix total 


25.6±2.5±3.2±6.1 



1.80±0.09±0.19±0.37±0.02 
2.18±0.23±0.33±1.15±0.03 
2.94±0.19±0.53±0.92±0.04 
0.86±0.10±0.16±0.31±0.01 
1.39±0.08±0.16±0.35±0.02 
1.62±0.21±0.41±1.21±0.02 
2.26±0.22±0.34±0.58±0.03 



2.98±0.19±0.53±0.93±0.04 
1.02±0.12±0.18±0.36±0.03 



TABLE IV: Systematic uncertainties on the signal yield and fit fractions. 



Efficiency qq 

DP PDF 



BB 
DP PDFs 



CR TOES & AE 
PDF parameters 



SCF rriES & AE 
PDF parameters 



BB niES & 
AE PDFs 



SCF qq niES 

fraction PDF parameters 



Fit bias Total 



Signal yield 


6.3 


19 


24 


35 


0.7 


22 


302 


2.1 


46 


310 


i?J (2460)^71+ FF 


0.0018 


0.0036 


0.0051 


0.00072 


0.00005 


0.0033 


0.010 


0.00009 


0.0037 


0.013 


D5 (2400)- 71+ FF 


0.003 


0.016 


0.024 


0.00063 


0.00013 


0.0065 


0.00096 


0.00011 


0.0023 


0.030 


piriOfD" FF 


0.016 


0.028 


0.031 


0.00069 


0.00010 


0.025 


0.0078 


0.00023 


0.0014 


0.052 


/2(1270)D° FF 


0.0054 


0.0091 


0.0077 


0.00040 


0.00005 


0.0078 


0.0029 


0.00006 


0.0011 


0.016 


D*{20W)-iT+ FF 


0.00097 


0.0028 


0.0045 


0.00034 


0.00004 


0.0027 


0.0091 


0.00007 


0.0052 


0.012 


Dn NR FF 


0.015 


0.023 


0.022 


0.00052 


0.00010 


0.020 


0.015 


0.00008 


0.0036 


0.043 


K matrix total FF 


0.0057 


0.014 


0.015 


0.00075 


0.00017 


0.012 


0.018 


0.00008 


0.010 


0.032 



These sources of systematic uncertainty are summarized in Table IIVI The total is obtained by combining all sources 
in quadrature. 

We consider additional systematic effects on the values of the branching fractions. These are uncertainties on the 
differences between the efficiencies of selection requirements on data and MC for tracking (1.0%), particle identification 
(4.0%), the neural network cut (3.2%), and the number of BB pairs (0.6%). Furthermore, where we have divided by 
a daughter branching fraction in order to isolate the B decay branching fraction, any uncertainty in the world average 
value used in the division also contributes systematic uncertainty. 

An additional source of uncertainty in Dalitz-plot analyses arises due to the composition of the Dalitz plot. We 
consider the following sources of model uncertainty: 

• Fixed parameters of contributing amplitudes: 

We vary the masses and widths of all resonances described by RBW shapes according to the uncertainties of 
the world average values [26] (with the exception of the Do(2400) mass, which we vary by ±100MeV/c^ to 
account for the discrepancy in the measured masses of charged and neutral isospin partners). We vary the a 
parameter of the Dn nonresonant contribution within its uncertainty. We change the radius parameter of the 
Blatt-Weisskopf factors from its nominal value of 4 GeV"^ to both 3 GeV"^ and 5 GeV"^. 

• Alternative parameterisations: 

We use the Gounaris-Sakurai lineshape 43| as an alternative description for the yo(770) resonance. We replace 



the TF+TT S-wave K-matrix term with contributions used in the analysis of i? — > D^tt+tf by Belle [22|, namely 
a (described as in Ref. ^), /o(980) (described by the Flatte distribution [H]) and /o(1370) (RBW). To address 
the possible discrepancy between the data and the fit result at low values of m-, we replace the Dtt nonresonant 
contribution with a functional form proposed for a putative "dabba" state [46|. We have also performed a fit 
in which the background from £'*-(2010)7r+ events escaping the veto is treated as a separate (seventh) BB 
background category, and include the deviation in the results as a source of model uncertainty. 

• Additional possible contributions: 
We repeat the fit adding states to the model: w(782), /9(145G), D(2600) (both as scalar and vector) and D(2760) 
(vector) ji^l. 



A summary is given in Table IVl The total model uncertainty is obtained by combining all sources in quadrature. 
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TABLE V: Dalitz-plot model uncertainties on the signal yield and fit fractions. Refer to the text for details of the model 
variations. 





Mass & 


Dtv NR 


BW barrier 


p(770)" GS 


Dtt S-wave 


TT^TT S-wave 






width 


a 


radius 


lineshape 


"dabba" 






Signal yield 


44 


6.4 


11 


1.1 


14 


67 




D5(2460)'7r+ FF 


0.028 


0.0027 


0.020 


0.00007 


0.0019 


0.0052 




DS(2400)'7r+ FF 


0.061 


0.031 


0.0098 


0.00066 


0.099 


0.043 




pi770)°D° FF 


0.045 


0.0056 


0.042 


0.0010 


0.00012 


0.034 




/2(1270)D° FF 


0.018 


0.00061 


0.0060 


0.00058 


0.0040 


0.014 




D;(2010)-7r+ FF 


0.018 


0.0028 


0.015 


0.00076 


0.025 


0.0097 




Dn NR FF 


0.10 


0.024 


0.021 


0.0060 




0.026 




K matrix total FF 


0.023 


0.0075 


0.010 


0.0034 


0.038 








Add a; (782) 


Add p(1450) 


Add D(2600) 
(scalar) 


Add D(2600) 
(vector) 


Add D(2760) 
(vector) 


7 BB Cat. 


Total 


Signal yield 


11 


53 


0.62 


13 


1.5 


8.4 


100 


D5(2460)-7r+ FF 


0.00019 


0.0060 


0.00026 


0.011 


0.00063 


0.00088 


0.037 


DS(2400)'7r+ FF 


0.0019 


0.00091 


0.0085 


0.011 


0.0090 


0.00034 


0.129 


p{770fD° FF 


0.015 


0.047 


0.00075 


0.050 


0.0096 


0.0015 


0.100 


/2(1270)L»° FF 


0.00048 


0.010 


0.00004 


0.021 


0.0021 


0.00050 


0.034 


D;(2010)~7r+ FF 


0.0015 


0.0095 


0.00012 


0.0040 


0.0016 


0.00088 


0.037 


Dtt NR FF 


0.0053 


0.045 


0.0025 


0.065 


0.0036 


0.0012 


0.136 


K matrix total FF 


0.0070 


0.034 


0.0019 


0.018 


0.011 


0.00073 


0.061 



DISCUSSION 

I sosp in symmetry can be used to relate the decay amplitudes of B^ 
SIS El: 



D°p°, B° -^ D-p+ and B+ -^ D°p+ [11 



(7) 
(8) 
(9) 

where A^/2 a-nd A1/2 are the amplitudes for isospin 3/2 and 1/2 final states respectively. These equations give the 
triangle relation 



A{D'p+) = 


\/3A3/2 , 




A{D-p+) = 


Vl73^3/2 + ^ 


/273A1/2 , 


v/2yl(^°p°) = 


V473A3/2 - ^ 


yv^A,/2 , 



A{D'^p+) = A{D-p+) + \/2A(D"p") ■ 



(10) 



This relation can be used to determine cos 6 Dp, where Sop is the phase between the ^3/2 and ^1/2 amplitudes, 
and Rop = |^i/2/>/2^3/2|- In QCD factorization, both of these are expected to be unity up to corrections due to 
final state interactions of O^Aqcb/itiq), where ^qcd is the QCD scale and mg is either rric or TUh [ll|. We 
constraints on these parameters using the same approach previously used in the D^*>tt system [SOJ, |5l|- Usi 



result for B{B'^ — )• £)°p°), together with world average values of B{B'^ 
lifetimes t{B+)/t{B°) ^, we find 



D-p+), B{B+ -> D"p+) and the ratio of 



cos Sdp 
Rdp 



o.998;°:j33 , 



where all sources of uncertainty are combined. These results suggest the presence of non-factorizable final state 
interaction effects that, in contrast to the D^*'7t system, do not introduce a significant non-zero phase difference 
between the isospin amplitudes. 

SUMMARY 

We have performed a Dalitz-plot analysis of 5° — ?► D^tt+tt" decays using the whole BaBAR dataset of 
470.9 ± 2.8 million BB events. We measure the inclusive branching fraction 



B{B° -^ D^TT+TT-) = (8.81 ± 0.18 ± 0.76 ± 0.78 ± 0.11) x 10" 
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where the first uncertainty is statistical, the second is systematic, the third is due to the Dalitz-plot model, and 
the fourth is due to secondary branching fractions. We find the Dalitz plot to be composed of contributions from 
1)2(2460)^, I?q(2400)~, p(770)° and /2(1270) as well as a tt+tt" S-wave, a Dn nonresonant S-wave term and a virtual 
£'*(2010)~ contribution. Wc determine their branching fractions: 

B{B^ -> i:>|(2460)-7r+) x B{Dl{2Am)- -^ Wt^-) = (1.80 ± 0.09 ± 0.19 ± 0.37 ± 0.02) x lO"'' , 

S(B° -^ Dl{2Am)-'r:+) x S(i:>*(2400)- -^ Wt^-) = (2.18 ± 0.23 ± 0.33 ± 1.15 ± 0.03) x lO"'' , 

B{B° ^ p{770fD°) = (2.98 ± 0.19 ± 0.53 ± 0.93 ± 0.04) X 10"\ 

B{B" ^ f2{1270)D^) = (1.02 ± 0.12 ± 0.18 ± 0.36 ± 0.03) X lO""*. 

Our Dalitz plot model differs from that obtained in a previous study of B^ -^ ZJ^tt+tt" by Belle [22| in that (i) 
we use the K-matrix description of the 7r"'"7r~ S-wave, instead of including separate contributions from the /o(600) 
(cr), /o(980) and /o(1370) scalar resonances; (ii) we include an additional Dtt nonresonant S-wave term. Our results 
for the inclusive branching fraction and for the color-suppressed decays i?° — > /3(770)°D° and B'^ — ?> /2 (1270)13'' 
are consistent with those from Belle and (for p(770)"Z3°) with theoretical predictions [13;|52[. However, we find the 
product branching fractions for the broad and narrow D** states (I3o(2400) and D2(2460), respectively) to have 
similar values. This result disagrees with the analysis by Belle, which found a much smaller value for the Dp (2400) 
branching fraction. 
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Appendix: K-matrix description of tt n S wave 

The K-matrix formalism gives a physical description of broad overlapping states - i. e. it does not violate unitarity, 
unlike the more conventional "sum of Breit-Wigners" approach. The K-matrix formalism can be shown to reduce to 
more familiar forms (the Breit-Wigner lineshape for single resonances, the Flatte lineshape [45| for coupled channels, 
the LASS formula 53] for broad resonances interferin g w ith nonresonant terms). Detailed descriptions of the K-matrix 



formalism can be found in various references 32h341 l54l[. Here we give an outline of the salient features. 

The scattering ("S") matrix describes transitions from initial states \i) into final states |/) {Sif = {f\S\i)), and can 
be written 

s = i+2t{p^y/^T{py^\ (A.ii) 

where / is the identity matrix, p is a diagonal phase-space matrix, with elements pu = 2qi/m with qi the threshold 
momentum, and T is the transition matrix. The unitarity requirement {SS' = S^ S = I) gives 

(r-i + ip)^ =r-i + zp, ^ K-^=T-^+ip, (A.12) 

where iiT is a Lorentz- invariant and Hermitian matrix which describes the decay process. This formalism was developed 
for scattering processes, but can also be applied to Dalitz plot analyses, with the assumption that the two "scattering" 
products do not interact with the third bachelor particle [33|. However, it is also necessary to include a process- 
dependent production vector, which accounts for the relative production rates of the different states \i). We refer to 
the "K-matrix amplitude" as a product of the production vector P and the (matrix) propagator (/ — iKp) : 

A,^{I-iKp)-lP, (A. 13) 



The K matrix is expressed as 
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K,,{s) 



„scatt 
■^0 



E 



(a) (a) 

9i 9) 



m-i. 



1 - SAQ 
S ~ SAO 



SA-ml 



(A.14) 



where the factor g^ is the real couphng constant of the K matrix pole a (with mass nia) to meson channel i, the 
parameters /f'^'^** and Sg'^'^** describe a smooth part for the K-matrix elements, and the last factor accounts for the 
so-called "Adler zero" , and suppresses kinematically fake singularities near tt+tt" production threshold (s represents 
the square of the tt+tt" invariant mass). The K-matrix parameters are determined from global fits to scattering data 
experiments below 1900MeV/c^ [3J]. Note that the phase-space for B^ — > D^ tt+tt^ extends beyond this limit, and 
that the K-matrix amplitude in this high-Tr+Tr" invariant mass region is therefore an extrapolation. 

The parameters unique to the production vector, by contrast, must be determined from our data. The P vector is 
given by 



P,(s) 



/'piod 



1 



prod 



prod 



E 



(3a9 



(q) 



(A.15) 



where as before the first term in the square brackets is nonresonant-like ("slowly varying"), and the second term 
is resonant-like. Hence the free parameters in the Dalitz plot fit are the complex coupling and production vector 
parameters Pa and /f™ (we use a fixed value of Sq"^" ). The index j runs over the open channels for the tttt S'-wave, 
which are: tttt, KK, rjrj^ rjrj' and 47r (or multi- meson) . At higher masses there are in principle more open channels, 
but this is not expected to affect the results significantly. Global fits to the scattering data determine the number 
of poles and their parameters. We use a 5 pole approximation, and give the values of all fixed parameters in the 
K-matrix model in Table IVll Note that all /P'° = for i 7^ 1 since we are interested only in the tttt final state. 



TABLE VI: K-matrix parameters from a global analysis of the available tttt scattering data from threshold up to 



1900 MeV/c^ 



Masses and coupling constants are given in GeV/c 



nia 




9kk 


a. 


g^r, 


a 


0.65100 


0.22889 


-0.55377 


0.00000 


-0.39899 


-0.34639 


1.20360 


0.94128 


0.55095 


0.00000 


0.39065 


0.31503 


1.55817 


0.36856 


0.23888 


0.55639 


0.18340 


0.18681 


1.21000 


0.33650 


0.40907 


0.85679 


0.19906 


-0.00984 


1.82206 


0.18171 


-0.17558 


-0.79658 


-0.00355 


0.22358 




rscatt 
Jll 


rscatt 
/l2 


/■scatt 
/l3 


rscatt 
/l4 


/■scatt 
/l5 




0.23399 


0.15044 


-0.20545 


0.32825 


0.35412 




^scatt 

■So 


prod 
^0 


Sao 


SA 






-3.92637 


-3.0 


-0.15 


1 





